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Abstract Single Degenerate model is the most widely accepted progenitor model of 
type la supernovae (SNe la), in which a carbon-oxygen white dwarf (CO WD) accretes 
hydrogen-rich material from a main sequence or a slightly evolved star (WD +MS) to 
increase its mass, and explodes when its mass approaches the Chandrasekhar mass limit. 
During the mass transfer phase between the two components, an optically thick wind 
may occur and the material lost as the wind may exist as circumstellar material (CSM). 
Searching the CSM around progenitor star is helpful to discriminate different progenitor 
models of SNe la. Meanwhile, the CSM is a source of color excess. The purpose of this 
paper is to study the color excess produced from the single-degenerate progenitor model 
with optically thick wind, and reproduce the distribution of color excesses of SNe la. 
Meng et al. (2009) systemically carried out binary evolution calculation of the WD +MS 
systems for various metallicities and showed the parameters of the systems before Roche 
lobe overflow and at the moment of supernova explosion in Meng & Yang (2009). With 
the results of Meng et al. (2009), we calculate the color excesses of SNe la at maximum 
light via a simple analytic method. We reproduces the distribution of color excesses of 
SNe la by our binary population synthesis approach if the velocity of the optically thick 
wind is taken to be of order of magnitude of 10 km s _1 . However, if the wind velocity is 
larger than 100 km s~\ the reproduction is bad. 
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1 INTRODUCTION 

Although type la supern ovae (SNe la) showed their importance i n determining cosmological param- 
eters, e.g. J7m and J7a dRiess et all 1 19981 iPerlmutter et all 1 19991). the progenitor system s of SNe la 
have not been confidently identified yet (IHillebrandt & Niemeverll2000tlLeibundgutll2000l) . It is widely 
believed that a SN la is from thermonuclear runaway of a carbon-oxygen white dwarf (CO WD) in 
a binary system. The CO WD accretes material from its companion to increase its mass. When its 
mass reaches its maximum stable mass, it explode s as a thermon uclear runaway and almost half of the 
WD mass is converted into radioactive nickel-56 dBranchll2004l) . Two basic scenarios have been dis- 
cussed over the last th ree decades. One is a single degenerate (SD) model, which is widely accepted 
(IWhelan&Ibenlll973l) . In this model, a CO WD increases its mass by accreting hydrogen- or helium- 
rich matter from its companion, and explodes when its mass approaches the Chandras ekhar mass limit. 
The companion may be a main-sequence star (WD+MS) or a red-giant star (WD+RG) dYungelson et al.l 
sunn & van den Heuvelll997l:lHachisu et al.ll999allHachisu et al.ll999bT:lNomoto et alll999Ll2003b 
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Lang er et al 1l2000l:lHan & Podsiadlowskil2004tlChen &D2007tlHanl2008HMeng et al.ll2009l:lLu et al 



2009) . An alternative is the theoretically less favored double degenerate (DD) model dlben & Tutukov 



1 984 IWebbinl 1984 . in which a system consisting of two CO WDs loses orbital angular momentum 



by gravitational wave radiation and merges finally. The merger m ay explode if the total mass of the 
system exceeds th e Chandrasekhar mass limit (see the reviews by iHillebrandt & Niemeverll2000l and 
iLeibundgutCOOOb . In theory, a large a mount of circumstel lar materials (CSM) may be produced via an 
optically thick wind for the SD model dHachisu et all 19961) . while there is no CSM around DD systems. 
Then, a basic method to distinguish the two progenitor models is to find the CSM around progenitor 

systems. 

Evidence for CSM was first found in SN2002ic dHamuy et alJ 120031) . which has shown ex- 
tremely pronounced hydrogen emission lines that have been interpreted as a sign of strong inter- 
action between supernova e jecta and CSM. The discovery of SN2002ic may uphold the SD model 
(lHan & Podsiadlowskil2006l). Recen t ly, the evidence for CSM was found in a normal SN la (SN 2006X) 
defined bv lBranch. Fisher & Nugent! (1 1993b and the CSM is proposed to be from a wind from a red-giant 



companion dPatat et al.ll2007al) . The C SM may pla y a key role to solve the problem of the low value 
of reddening ratio of external galaxy dWang||2005l) . which is very important for precision cosmology 
dWang et al.ll2008h . 

If a SN la is surrounded by a large amount of CSM, its color observed shou ld be redder than 
its intrinsic color, which results in a color excess, E(B - V). iReindl et al.l d2005l) showed the color 
excesses of more than one hundred SNe la at maximum light, which suggests a mission to check which 
progen itor model of SNe la can explain the distribution of the color excesses. Recently, iMeng et ail 
d2009l) performed binary stellar evolution calculations for more than 25,000 close WD binary systems 
with various metallicities, and present all the parameters of the systems for SNe la be fore the Roche 
lobe o verflow (RLOF) and at the mome nt of supernova explo sion in a following paper dMeng & Yand 
120091) . In their works, the prescription of lHachisu et al.ldl999al) for the accretion efficien cy of hydrogen- 
rich material was adopted by assuming an optically thick wind dHachisu et al.l 1 1996b . and then their 
works provide a possibility to check whether the SD model with optically thick wind can reproduce the 
distribution of color excesses of SNe la obtained from observatio n or not. The purpo se of this paper is 
to check the possibility, and this work is based on the results from lMeng et al.l d2009l) . 

In section [2] we describe our model. We show the results in section [3] and give discussions and 
conclusions in sections |4] and [5] 



2 MODEL AND PHYSICS INPUTS 



2.1 the mass distribution of lost hydrogen-rich material 

As described in section Q] the color of a supernova is reddened by the CSM surrounding the supernova. 
We first check whether there exists enough CSM resulting from op tically thick wind. The following is 
a simple description about how to obtain the CSM. As described in lMeng et al.l d2009l) . in a WD + MS 
system, the companion fills its Roche lobe at MS or during HG, and transfers material onto the WD. 
If the mass-transfer rate, M2I, exceeds a critical value, M cr , we assume that the accreted hydrogen 
steadily burns on the surface of WD and that the hydrogen-rich material is converted into helium at the 
rate of M cr . The unprocessed matter is assumed to be lost from the system as an o ptically thick wind at a 
rate of M w ind = |M 2 | - M cr dHachisu et al.ll 1996b . Adopting the prescription oflHachisu et al.ldl999ab 
on WDs accreting hydrogen-rich material from their companions. IMeng et al.l d2009h obtained the initial 
parameters of WD + MS systems for SNe la, and at the same time, the final parameters at the moment 
of SN la explosion, such as Mf N . In this paper, incorporating the binary evolution results in lMeng et al.l 
d2009b into the rapid binary evolution code developed by Hurley et al. (2000, 2002), we carry out a series 
of binary population synthesis (BPS) studies for various Z. In each BPS study, 10 7 binaries are generated 
by Monte Carlo simulation and a circular orbi t is ass umed for all binaries. The basic parameters for the 
simulations are same to those in IMeng et al l ( l2009h . It can be shown that a WD+MS system may 
originate from three possible evolution channels, namely, the He star channel, the EAGB channel 
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and the TPAGB channel (see iMeng et alj|2009l for details). We assume that a SN la is produced 
if the initial parameters of a WD + MS system, e.g. initial orbital period P^. h and initial secondary 
mass M\, locate in the appropriate regions of the parameters for SNe la at the onset of RLOF. We can 
obtain the companion mass at the moment of explosion by interpolation in the thr ee-dimensional gri d 
(M^ D , M^logP 1 ) of the more than 25,000 close WD binary system calculated in lMeng etall d2009l) . 
In Fig. Q] we show the mass distribution of hydrogen-rich material lost as the optically thick wind for 
various metallicities, M W j n( j = [M\ + ^wd) — C^f N + ^wd)' where superscript i and SN represent 
the initial and final values for white dwarf and secondary, respectively, and is assumed to be 

1.378M Q . We can see from the figure that the distribution of the lost mass peaks at about 0.3 Mq and 
has a high-mass tail. The amount of the lost material may be as large as 2.5 Mq, which should contribute 
to the color excess of SNe la. 

In Fig.Q] we can see that there does not seems to be a systemic trend with metallicity. Actually, the 
influence of metallicity on the M W j n d is complicated. M w i n d is mainly determined by M\ and A£\y D . 
The two parameters are both systemically affected by metallicity, but the tendency is reversed, i.e. the 
peak of the distribution of companion m ove to higher mas s with metallicity, while the peak for WD 
mass to lower mass (see Figs. 9 and 10 in lMeng et al.ll2009b . In a ddition, the metallic ity also affect the 
mass-transfer rate between WD and its companion, and then Mf N (lLanger et al.l2000T) . The complicated 
influence of metallicity on the distribution of A/ W i n d results in a non- systemic trend of the distribution 
with metallicity. We also noticed that the percentage of high M w ; n d, e.g. M W j n d > 2.0M Q , increases 
with metalliticity. High M w i n d is mainly determined by the upper boundary of the companion mass, 
which moves to higher mass with metallicity (see Fig. 4 in IMeng et al.l l2009h . So, a binary system 
producing SN la with high metallicity may loss more hydrogen-rich material by optically thick wind, 
and then the high-metallicity model in Fig. [T]shows a higher percentage of high Af w ; n d. 
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Fig. 2 The distribution of color excesses of SNe la at maximum light for various metallicties, 
where wind velocity is assumed to be 1 km s -1 . The solid histogram is from observation, and 
the bar represents its maximum error dReindl et al.ll2005h . The black dashed line is the sum 
of those of Z=0.01, 0.02 and 0.03, where the weights for the three components are 20%, 40% 
and 40% respectively. Via K-S test, the percentages indicated in the figure show the confident 
level that the distributions of color excess from theory and observation are indistinguishable. 



2.2 model 

iBohlin et all dl978l) found 

E(B - V)/n K = 1.72 x 10" 22 mag cm 2 , (1) 

for th e Galaxy, whe re nn is total hydrogen column density and this relation linearly depends on metal- 
licity dDrainell2003l) . The linear relation can be fitted by 

E(B - V)/n H = (17.4 x Z/Z Q - 0.454) x 10~ 23 mag cm 2 , (2) 

where the fitting data are from lBohlin et al.ldl978b . lKoorneefl(ll982l) . lFitzpatrickl(ll985b and lMartin et al.l 
In the following, we describe how to obtain the total hydrogen column density, nu- 
While the CSM around the progenitor of a SN la may be asymmetric, for simplicity, we assume 
that the mass loss of optically thick wind is spherically symmetric and the mass-loss rate is constant 
during the whole mass-transfer phase. Then, hydrogen number density is only a function of the distance 
of wind material from the progenitor star, i.e. 

n(r) = ar~ 2 , (3) 

where a is a coefficient determined by 

f ar~ 2 ■ 4nr 2 dr = 7V H , (4) 
Jr 
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where rrj is the radius of progenitor, r m is the maximum distance which hydrogen-rich material can 
arrive as optically thick wind at the moment of SN la explosion, and Nn is the total number of hydrogen 
atom. The total hydrogen column density at maximum light is calculated from 

ar~ 2 dr = tih, (5) 



where r x is the distance of explosion ejecta from explosion center at maximum light. From equations 
(4) and (5), we can obtain 

n H = ■ (6) 



A^h can be obtained from 



iV H = I™!*!*., (7) 

TOH 



where M w i n( j is the total mass lost as optically thick wind, mn is the mass of a hydrogen atom and X is 
the mass fraction of hydrogen and is set to be 

X = 0.76 - 3.0Z, (8) 

dPols etalJll998h . r m can be obtained from 

^wind^d; (9) 

where Vynnd is the velocity of the optically thick wind, and is the delayed time from the onset of 
mass transfer to the moment of SN la explosion. From the results of IMeng et al.l d2009l) . td can be 
approximated by 

log(t d /yr) = -j?Mj + 7.8, (10) 

where is the initial mass of the secondary (the mass donor unit in solar mass) in a WD + MS system. 
td from the equation is a rough estimation of the mean value for a certain M\ and has an error of about 
50%. n, the distance of explosion ejecta from explosion center at maximum light, can be obtained by 
the product of the velocity of supernova ejecta an d the rise time of lig ht curve of SNe la. We simply 
assume that the velo city of ejecta is 10000 km s" 1 ( Gamezo et al.ll2003l) and the rise time of light curve 



of SNe la is 20 day dConlev et alj|200d: IStrovinkll2007t) . The ejecta velocity adopted here corresponds 
to a typical p hotospheric velocity dWang et al.ll2.Q03r) and might be lower than terminal ejecta velocity 
dWandl2006l) . However, the uncertainty resulting from the ejecta velocity is moderate and accepted. We 
will discuss its influence in section[4] The rise time does not significantly affect the final results. 

IMeng et all d2009l) performed binary stellar evolution calcu l ations for more than 25,000 close 
WD binary systems with various metallicities and IMeng & Yand (|2009|) presented the distribution of 
all the parameters for these close systems before the RLOF and at th e moment of SN la explo sion. 
Incorporating their results into the binary population synthesis code of iHurlev et al.l d2000l |2002|) . we 
obtain the distribution of the wind mass, M w j„d (see Fig.Q]) and the c olor excess via equ ations[T]and|6l 
The basic parameters for Monte Carlo simulations are same to that in lMeng et alj d2009l) when primor- 
dial binary samples are generated. Because the code is valid just for Z < 0.03, only five metallicities 
(i.e. Z = 0.03, 0.02, 0.01, 0.004 and 0.001) are examined here. 

The greatest uncertainty of our model is from Vwind- Here, we assume that Kvi n d = 10 km s _1 and 
we will discuss it in section |4] 
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Fig. 3 The distribution of color excesses of SNe la at maximum light for Z = 0.02 and 
various wind velocities (in km s" 1 ). The solid hi stogram is from observation, and the bar 
represents its maximum error dReindl et al.l 12005*1) . Via K-S test, the percentages show the 
confident level that the distributions of color excess from theory and observation are indistin- 
guishable. The black dashed line is the best fitted line, where the weights for the components 
of Z=0.01, 0.02 and 0.03 are 20%, 10% and 70%, respectively. 



3 RESULTS 

In Fig. |2] we show the distribution of the color excesses of SNe la at maximum light for v arious metal- 
licitie s. The distribution from observation is also shown by solid histogram in the figure dReindl et al.l 
120051) . We see in the figure that the SD model with optically thick wind can reproduce the distribution 
of color excesses of SNe la observed. The K-S test shows that the cases of Z = 0.02 and Z = 0.03 
have the highe st confident level th at the distributions from theory and observation are indistinguishable. 
The sample of lReindl et alj d2005l) includes the SNe la with various metallicities. H owever, it is difficul t 
to determine the m etallicity of the host galaxy of every SN la for a large distance dHamuv et al.ll2000l) . 
IWang et ail d2006l) collected the properties of several SNe la from previous papers and noticed that the 
mean metallicity of the host galaxies of the SNe la is [12 + log(0/H)] moan = 8.85 ±0.10 (private 
comm unication), which is consistent with solar metallicity ([12 + log(0/H)] so i ar ~ 8.8 dZaritskv et al 



1994). Then, our result is consistent with observation although the sample collected by I Wang et al 
d2006h is small. 

There are some SNe la w hose color e xcesses are very large (i .e. larger than 1. 0), such as SN199 9cl 
Jjhall2002t [Reindl et alj|2005l). SN 2003cg (lElias-Rosa et al.ll2006l) and SN2006X dWang et al.ll2008b . In 



the sample of Reindl et al.l ( 20051) . only one among 113 SNe la has a color excess larger than 1. In our 



simulation, the possibility of the high-color-excess SNe la is 2.0% for Z = 0.03, 0.4% for Z = 0.02 
and for Z < 0.02, consistent with that of the observations. This result might imply that the SNe la 
may have a color excess larger than 1 only when their host galaxies have a metallicity larger than 0.02. 
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Observationally, t he host galaxies of SN1999cl and SN2006X , i.e. NGC4501 and NGC4321, are both 
oversolar galaxy dCaputo et al.ll2000t iDors Jr & Copettill200r3) . Although there is no information about 
the metallicity of the h ost galaxy of SN200 3cg (NGC3169), NGC3169 is a Sa galaxy and probable has 
a oversolar metallicity dWillner et al.ll 19851) . 



4 DISCUSSION 



Our analytic model is so simple. We discuss various uncertainties about our model in this section. 



4.1 wind velocity 

The major uncertainty of our model is from the assumption that V w { n( i = 10 km s _1 . Although the 
consistency between the theoretical distribution of color excess and that from observation upholds 
this assumption, there i s no direct observational evidence to verify it. Many observational efforts were 
projected to find C SM dHamuv et al.ll2003l lAldering et aT1l2006t iPanagia et al1l2006t lOfek et alj|2007t 



Patat et al1l2.007 allbb . and only one observation obtained the constraint of the wind velocity ( Patat et aTl 
2007ah . The wind velocity is constrained to be smaller than 50 km s 1 (see Fig. 2 in lPatat et alj|20073h 



and then Patat et al. suggested that the progenitor of SN2006X should be a WD+RG system. However, 
we can not rule out the possibility that the observed CSM is fro m the opti cally thick wind since the 
companion has not been directly observed dHachisu et al. 2008). Recently, iBadenes et al.l d2007l) ex- 
plored the relationship between the SD models with optically thick wind for SNe la and the properties 
of the supernova remnants that evolve after the explosion. They found that the optically thick wind with 
velocity larger than 200 km s" 1 would excavate large low-density cavities around the progenitors. The 
large cavities are incompatible with the dynamics of the forward shock and the X-ray emission from the 
shocked ejecta in all the SNe la remnants as they examined. However, they also showed that if a wind 
velocity of 10 km s _1 is adopted, the properties of type la supernova remnants are well compatible with 
the prediction from the SD model with optically thick wind. Generally, the escape velocity from a white 
dwarf is in an order of m agnitude of 10 3 km s _1 , which is upheld by observations from rec urrent novae 
dWood & Locklevl2000l) . In the theoretically framework laid down bv lHachisu et alj dl996|). an optically 
thick w ind is formed in CO WD envelope with photospheric velocity of ~ 10 3 km s" 1 dHachisu et al.l 
Il999albh . where the expansion o f the photosphere is driven by helium flash in a helium shell on top of 
CO WD dKato&Hachisulll999t) . Then, in both theory and observation, it seems not to uphold an as- 
sumption of a low wind velocity at present. We check the influence of wind velocity on the distribution 
of color excess for the case of Z = 0.02 and the results are shown in Fig. [3] We see in the figure that 
when V > 100 km s~\ K-S test gives a low confident level that the distributions of color excess from 
theory and observation are indisting uishable, while when V < 50 km s _1 , K-S test shows an acceptable 
level. This result is similar to that in lBadenes et al.l d2007l) . 



4.2 delay time 

The secondary uncertainty of our model is from t^. Since optically thick wind may stop before SN la 
explosion, used in this paper overestimates the delayed time of the wind for some systems. However, 
this is not a serious pro blem since for most cases, SNe la occur during t he wind phase or after the wind 
phase for a short time dHan & Podsiadlowskil 120041; iMeng et ai1l2009l) . In addition, obtained from 
equation (10) only approximates the mean value for a given M\. This is also not a serious problem since 
we only check an average distribution of color excess. Generally, the peak of the distribution of color 
excesses moves to lower color excess with t^, and the peak value increases while the percentage of high 
color-excess SNe la decreases with t^. 
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4.3 metallicity 

The combination of different metallicities is checked and shown in Fig.[2]by the black dashed line since 
not all SNe la are from a unique metal licity. The we i ghts o f the components for Z= 0.01, 0.02 and 0.03 
are simply fitted from the samples of lHamuy et al.l ( 2000l) and IWang et al.l (120061) . We also show the 
best fitted line in Fig.[3]by the black dashed line. We can see from these lines that the combination of 
different metallicities does not significantly increase the confident level that the distributions for theory 
and observation are indistinguishable. Actually, for any combination of Z = 0.01, 0.02 and 0.03, the 
results will be acceptable, i.e. yielding similar K-S test probability, which is derived from similar K-S 
test probability for the distributions of Z = 0.01, 0.02 and 0.03. So, the basic results here still hold. 

4.4 interstellar extinction 

While we are not su re about the existen ce of the optically thick wind (it is a prediction by a model, 
not an observed fact, lHachisu et aill 19961) . we are absolutely sure about the existence of the interstellar 
extinction (IE). In the paper, we did not consider the influence of IE since it is difficult to separate it from 
CSM dust. However, we know that Type la SNe, as any other celestial object, suffer from interstellar 
reddening, arising in material that has nothing to do with the circumstellar environment of the exploding 
star. Let's consider the case of SN 2006X, whose color excess is certainly larger than E(B-V) =1, while is 
most likely arising in a cold molecular cloud which has nothing to do with the explosion site dWang et al.l 
120081) . So, for the case suffering interstellar reddening, we should remove the influence of IE on color 
excess. But at present, it is very difficult to do this. However, it is a clear effect that IE will lead to certain 
directions of where parameters such as Vwind and V c j cc should go, e.g. moving some very reddened 
observed SNe la to lowering reddened one would yield a higher Vvvind or y c j CC . For example, if we 
assume rather arbitrarily and simply that all observed SNe la suffer a extinction of E(B — V^)host = 0.1 
within its host galaxy, our model suggests Vwind = 20 km s _1 or = 20000 km s . Similarly, if 
E(B — y)host = 0.2, our model suggests Vwind = 50 km s _1 or V^j ec = 50000 km s _1 . So, the effect 
of IE could be counteracted by another uncertainties, V e j e c, an d tne basic results in this paper are still 
valid, (see the discussion about ejecta velocity below). 

If the color excess is from CSM as shown in this paper, one may argued that we would have seen 
the CSM material through radio and X-ray emission, ar ising in the shock produced by the fast moving 
SN ejecta crashing into the CSM dStockdale et alj[2~006l) . However. iPanagia et al.l (|2006j) showed a non- 
detection result for radio for 27 SNe la via the Very Large Array (VLA) observations, even including 
SN 2002ic, which indicated that mass-lose rate should be lower than ~ 3 x 10~ 8 M^yr . We take 
a typical case to chec k whether our mo del contradicts with observations. We set companion mass to 
2.1M Q (see Fig .10 in lMeng et al]|2009l) . which corresponds to a delay time of 2.5 x 10 6 yr. The mass 
of lost material is set to O.3M (see Fig.Q]). The mean mass-lose rate is 12 x 10~ 8 M Q yr _1 , which is 
higher but still comparable to that inferred fr om observations. A t tempts to detect radio emission from 
SN 2002ic with the VLA were unsuccessful dBerger et alj|2003t IStockdale etafll2003l) . SN 2002ic is 



the first case to show the signal of CSM flH amuy et al.ll2003l) and the amount of the CSM is 0.5-6M, 



dWang et al.ll200llChugai & Yungelsonll2004HUenishi et al.ll2004t iKotak et ai]|2004j) and its properties 



may be explained successful by the WD+MS mo del used here ( Han & Podsiadlowski 20061) . In addition, 
another r ecent twins of SN 2002ic (SN 2005gj. lA.ldering et aL 20061) has also not been detected using 



the VLA dSoderberg & Fraill2005i) . These non-detecti on results may indic ate that the mechanism that is 
successful used for SN Ib/c may not work for SNe la dPanagia et al.ll2006h . So, our model is at least not 
inconsistent with observation at present. 

4.5 ionization and ejecta velocity 

Generally, when a SN explodes, its radiation field is quite strong, especially at distance rjj < 10 16 
cm (for SN 2006X, rjj < 4 x 10 15 , |Patat et alj|2007al) . This radiation can easily ionize hydrogen and 
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Fig. 4 The distribution of color excesses of SNe la at maximum light for different ejecta 
velocity, where wind velocity is assumed to be 10 km s _1 . The ejecta velocity is unit in 10 4 
km s _1 . The black dashed line is the best fitted line, where the weights for the components of 
Z=0.01, 0.02 and 0.03 are 20%, 10% and 70%, respectively. 



evaporate dust up to quite large distances, causing the disappearance of reddening. In our simple model, 
this effect is not considered completely. We check the influence of the effect. 

We set ?'h = 10 16 cm and assume that all hydrogen atoms are ionized within the shell with radius 
of rn- The column density of ionized hydrogen is ?ih,o = / ri H ar~ 2 dr, where r\ = 1.73 x 10 15 cm (a 
ejecta velocity of 10 4 km s _1 and a rise time of 20 day are assumed). Then, the relative uncertainty from 
the ionization effect is = ^ H ' ~ r , l)rm ~ rH '°~ ri ~ 0.83 (r m > n). The uncertainty seems too 

>ih (r m -ri)r H ,o m,o 

large. However, this situation may be improved by taking a higher ejecta velocity and a smaller distance 
rn- A distance of rjj = 4 x 10 15 like SN 2006X may reduce the uncertainty to 0.57. 

The ejecta velocity used here ma y be smaller than the terminal velocity of SN la which can be as 
high as 3 x ^kms" 1 dWandl2006l) . This value may reduce the uncertainty from the ionization effect 
to ~ 0.48, which means that no more than a half of CSM are evaporated up to large distances and have 
no influence on the color excess of SNe la. Then, the ionization effect at most decreases the confident 
level that theory and observation are indistinguishable from 74% to 65%. 

In Fig. |U we show the influence of ejecta velocity on the distribution of the color excess of SNe 
la. The influence is similar to that of wind velocity. From the figure, we can see that the peak of the 
distribution increases with the ejecta velocity, and when V c j cc < 5 x 10 4 km s _1 , the distributions of 
color excess from theory and observation are indistinguishable on the confident level of higher than 
50%. Fortunately, the effect of ionization and a higher ejecta velocity is same to that of interstellar, and 
their influence on the decrease of the confident level may counteract each other. Then, we can say that 
the distribution of color excess between theory and observation are indistinguishable on the level of 
higher than 59%, even when y c j CC = 3 x 10 4 km s _1 . 
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For the discussion above, our model may get a meaning result although our analytic model is so 
simple. 

5 SUMMARY AND CONCLUSIONS 

In summary, if a wind velocity of ~ 10 km s _1 is adopted, the SD model with optically thick wind may 
reproduce the distribution of color excess of SNe la obtained from observation, which might support the 
SD model for SNe la. However, if a wi nd velocity larger than 100 km s _1 is adopted, the reproduction is 
bad. Our results are similar to those of lBadenes et al.l d2007l) . Then, it should be encouraged to perform 
more detailed observations about the velocity of CSM around progenitors of SNe la. 
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Abstract Single Degenerate model is the most widely accepted progenitor model of 
type la supernovae (SNe la), in which a carbon-oxygen white dwarf (CO WD) accretes 
hydrogen-rich material from a main sequence or a slightly evolved star (WD +MS) to 
increase its mass, and explodes when its mass approaches the Chandrasekhar mass limit. 
During the mass transfer phase between the two components, an optically thick wind 
may occur and the material lost as the wind may exist as circumstellar material (CSM). 
Searching the CSM around progenitor star is helpful to discriminate different progenitor 
models of SNe la. Meanwhile, the CSM is a source of color excess. The purpose of this 
paper is to study the color excess produced from the single-degenerate progenitor model 
with optically thick wind, and reproduce the distribution of color excesses of SNe la. 
Meng et al. (2009) systemically carried out binary evolution calculation of the WD +MS 
systems for various metallicities and showed the parameters of the systems before Roche 
lobe overflow and at the moment of supernova explosion in Meng & Yang (2009). With 
the results of Meng et al. (2009), we calculate the color excesses of SNe la at maximum 
light via a simple analytic method. We reproduces the distribution of color excesses of 
SNe la by our binary population synthesis approach if the velocity of the optically thick 
wind is taken to be of order of magnitude of 10 km s _1 . However, if the wind velocity is 
larger than 100 km s~\ the reproduction is bad. 

Key words: Stars: white dwarfs - stars: supernova: general 



1 INTRODUCTION 

Although type la supern ovae (SNe la) showed their importance i n determining cosmological param- 
eters, e.g. J7m and J7a dRiess et all 1 19981 iPerlmutter et all 1 19991). the progenitor system s of SNe la 
have not been confidently identified yet (IHillebrandt & Niemeverll2000tlLeibundgutll2000l) . It is widely 
believed that a SN la is from thermonuclear runaway of a carbon-oxygen white dwarf (CO WD) in 
a binary system. The CO WD accretes material from its companion to increase its mass. When its 
mass reaches its maximum stable mass, it explode s as a thermon uclear runaway and almost half of the 
WD mass is converted into radioactive nickel-56 dBranchll2004l) . Two basic scenarios have been dis- 
cussed over the last th ree decades. One is a single degenerate (SD) model, which is widely accepted 
(IWhelan&Ibenlll973l) . In this model, a CO WD increases its mass by accreting hydrogen- or helium- 
rich matter from its companion, and explodes when its mass approaches the Chandras ekhar mass limit. 
The companion may be a main-sequence star (WD+MS) or a red-giant star (WD+RG) dYungelson et al.l 
sunn & van den Heuvelll997l:lHachisu et al.ll999allHachisu et al.ll999bT:lNomoto et alll999Ll2003b 
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Lang er et al 1l2000l:lHan & Podsiadlowskil2004tlChen &D2007tlHanl2008HMeng et al.ll2009l:lLu et al 



2009) . An alternative is the theoretically less favored double degenerate (DD) model dlben & Tutukov 



1 984 IWebbinl 1984 . in which a system consisting of two CO WDs loses orbital angular momentum 



by gravitational wave radiation and merges finally. The merger m ay explode if the total mass of the 
system exceeds th e Chandrasekhar mass limit (see the reviews by iHillebrandt & Niemeverll2000l and 
iLeibundgutCOOOb . In theory, a large a mount of circumstel lar materials (CSM) may be produced via an 
optically thick wind for the SD model dHachisu et all 19961) . while there is no CSM around DD systems. 
Then, a basic method to distinguish the two progenitor models is to find the CSM around progenitor 

systems. 

Evidence for CSM was first found in SN2002ic dHamuy et alJ 120031) . which has shown ex- 
tremely pronounced hydrogen emission lines that have been interpreted as a sign of strong inter- 
action between supernova e jecta and CSM. The discovery of SN2002ic may uphold the SD model 
(lHan & Podsiadlowskil2006l). Recen t ly, the evidence for CSM was found in a normal SN la (SN 2006X) 
defined bv lBranch. Fisher & Nugent! (1 1993b and the CSM is proposed to be from a wind from a red-giant 



companion dPatat et al.ll2007al) . The C SM may pla y a key role to solve the problem of the low value 
of reddening ratio of external galaxy dWang||2005l) . which is very important for precision cosmology 
dWang et al.ll2008h . 

If a SN la is surrounded by a large amount of CSM, its color observed shou ld be redder than 
its intrinsic color, which results in a color excess, E(B - V). iReindl et al.l d2005l) showed the color 
excesses of more than one hundred SNe la at maximum light, which suggests a mission to check which 
progen itor model of SNe la can explain the distribution of the color excesses. Recently, iMeng et ail 
d2009l) performed binary stellar evolution calculations for more than 25,000 close WD binary systems 
with various metallicities, and present all the parameters of the systems for SNe la be fore the Roche 
lobe o verflow (RLOF) and at the mome nt of supernova explo sion in a following paper dMeng & Yand 
120091) . In their works, the prescription of lHachisu et al.ldl999al) for the accretion efficien cy of hydrogen- 
rich material was adopted by assuming an optically thick wind dHachisu et al.l 1 1996b . and then their 
works provide a possibility to check whether the SD model with optically thick wind can reproduce the 
distribution of color excesses of SNe la obtained from observatio n or not. The purpo se of this paper is 
to check the possibility, and this work is based on the results from lMeng et al.l d2009l) . 

In section [2] we describe our model. We show the results in section [3] and give discussions and 
conclusions in sections |4] and [5] 



2 MODEL AND PHYSICS INPUTS 



2.1 the mass distribution of lost hydrogen-rich material 

As described in section Q] the color of a supernova is reddened by the CSM surrounding the supernova. 
We first check whether there exists enough CSM resulting from op tically thick wind. The following is 
a simple description about how to obtain the CSM. As described in lMeng et al.l d2009l) . in a WD + MS 
system, the companion fills its Roche lobe at MS or during HG, and transfers material onto the WD. 
If the mass-transfer rate, M2I, exceeds a critical value, M cr , we assume that the accreted hydrogen 
steadily burns on the surface of WD and that the hydrogen-rich material is converted into helium at the 
rate of M cr . The unprocessed matter is assumed to be lost from the system as an o ptically thick wind at a 
rate of M w ind = |M 2 | - M cr dHachisu et al.ll 1996b . Adopting the prescription oflHachisu et al.ldl999ab 
on WDs accreting hydrogen-rich material from their companions. IMeng et al.l d2009h obtained the initial 
parameters of WD + MS systems for SNe la, and at the same time, the final parameters at the moment 
of SN la explosion, such as Mf N . In this paper, incorporating the binary evolution results in lMeng et al.l 
d2009b into the rapid binary evolution code developed by Hurley et al. (2000, 2002), we carry out a series 
of binary population synthesis (BPS) studies for various Z. In each BPS study, 10 7 binaries are generated 
by Monte Carlo simulation and a circular orbi t is ass umed for all binaries. The basic parameters for the 
simulations are same to those in IMeng et al l ( l2009h . It can be shown that a WD+MS system may 
originate from three possible evolution channels, namely, the He star channel, the EAGB channel 
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and the TPAGB channel (see iMeng et alj|2009l for details). We assume that a SN la is produced 
if the initial parameters of a WD + MS system, e.g. initial orbital period P^. h and initial secondary 
mass M\, locate in the appropriate regions of the parameters for SNe la at the onset of RLOF. We can 
obtain the companion mass at the moment of explosion by interpolation in the thr ee-dimensional gri d 
(M^ D , M^logP 1 ) of the more than 25,000 close WD binary system calculated in lMeng etall d2009l) . 
In Fig. Q] we show the mass distribution of hydrogen-rich material lost as the optically thick wind for 
various metallicities, M W j n( j = [M\ + ^wd) — C^f N + ^wd)' where superscript i and SN represent 
the initial and final values for white dwarf and secondary, respectively, and is assumed to be 

1.378M Q . We can see from the figure that the distribution of the lost mass peaks at about 0.3 Mq and 
has a high-mass tail. The amount of the lost material may be as large as 2.5 Mq, which should contribute 
to the color excess of SNe la. 

In Fig.Q] we can see that there does not seems to be a systemic trend with metallicity. Actually, the 
influence of metallicity on the M W j n d is complicated. M w i n d is mainly determined by M\ and A£\y D . 
The two parameters are both systemically affected by metallicity, but the tendency is reversed, i.e. the 
peak of the distribution of companion m ove to higher mas s with metallicity, while the peak for WD 
mass to lower mass (see Figs. 9 and 10 in lMeng et al.ll2009b . In a ddition, the metallic ity also affect the 
mass-transfer rate between WD and its companion, and then Mf N (lLanger et al.l2000T) . The complicated 
influence of metallicity on the distribution of A/ W i n d results in a non- systemic trend of the distribution 
with metallicity. We also noticed that the percentage of high M w ; n d, e.g. M W j n d > 2.0M Q , increases 
with metalliticity. High M w i n d is mainly determined by the upper boundary of the companion mass, 
which moves to higher mass with metallicity (see Fig. 4 in IMeng et al.l l2009h . So, a binary system 
producing SN la with high metallicity may loss more hydrogen-rich material by optically thick wind, 
and then the high-metallicity model in Fig. [T]shows a higher percentage of high Af w ; n d. 
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Fig. 2 The distribution of color excesses of SNe la at maximum light for various metallicties, 
where wind velocity is assumed to be 1 km s -1 . The solid histogram is from observation, and 
the bar represents its maximum error dReindl et al.ll2005h . The black dashed line is the sum 
of those of Z=0.01, 0.02 and 0.03, where the weights for the three components are 20%, 40% 
and 40% respectively. Via K-S test, the percentages indicated in the figure show the confident 
level that the distributions of color excess from theory and observation are indistinguishable. 



2.2 model 

iBohlin et all dl978l) found 

E(B - V)/n K = 1.72 x 10" 22 mag cm 2 , (1) 

for th e Galaxy, whe re nn is total hydrogen column density and this relation linearly depends on metal- 
licity dDrainell2003l) . The linear relation can be fitted by 

E(B - V)/n H = (17.4 x Z/Z Q - 0.454) x 10~ 23 mag cm 2 , (2) 

where the fitting data are from lBohlin et al.ldl978b . lKoorneefl(ll982l) . lFitzpatrickl(ll985b and lMartin et al.l 
In the following, we describe how to obtain the total hydrogen column density, nu- 
While the CSM around the progenitor of a SN la may be asymmetric, for simplicity, we assume 
that the mass loss of optically thick wind is spherically symmetric and the mass-loss rate is constant 
during the whole mass-transfer phase. Then, hydrogen number density is only a function of the distance 
of wind material from the progenitor star, i.e. 

n(r) = ar~ 2 , (3) 

where a is a coefficient determined by 

f ar~ 2 ■ 4nr 2 dr = 7V H , (4) 
Jr 
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where rrj is the radius of progenitor, r m is the maximum distance which hydrogen-rich material can 
arrive as optically thick wind at the moment of SN la explosion, and Nn is the total number of hydrogen 
atom. The total hydrogen column density at maximum light is calculated from 

ar~ 2 dr = tih, (5) 



where r x is the distance of explosion ejecta from explosion center at maximum light. From equations 
(4) and (5), we can obtain 

n H = ■ (6) 



A^h can be obtained from 



iV H = I™!*!*., (7) 

TOH 



where M w i n( j is the total mass lost as optically thick wind, mn is the mass of a hydrogen atom and X is 
the mass fraction of hydrogen and is set to be 

X = 0.76 - 3.0Z, (8) 

dPols etalJll998h . r m can be obtained from 

^wind^d; (9) 

where Vynnd is the velocity of the optically thick wind, and is the delayed time from the onset of 
mass transfer to the moment of SN la explosion. From the results of IMeng et al.l d2009l) . td can be 
approximated by 

log(t d /yr) = -j?Mj + 7.8, (10) 

where is the initial mass of the secondary (the mass donor unit in solar mass) in a WD + MS system. 
td from the equation is a rough estimation of the mean value for a certain M\ and has an error of about 
50%. n, the distance of explosion ejecta from explosion center at maximum light, can be obtained by 
the product of the velocity of supernova ejecta an d the rise time of lig ht curve of SNe la. We simply 
assume that the velo city of ejecta is 10000 km s" 1 ( Gamezo et al.ll2003l) and the rise time of light curve 



of SNe la is 20 day dConlev et alj|200d: IStrovinkll2007t) . The ejecta velocity adopted here corresponds 
to a typical p hotospheric velocity dWang et al.ll2.Q03r) and might be lower than terminal ejecta velocity 
dWandl2006l) . However, the uncertainty resulting from the ejecta velocity is moderate and accepted. We 
will discuss its influence in section[4] The rise time does not significantly affect the final results. 

IMeng et all d2009l) performed binary stellar evolution calcu l ations for more than 25,000 close 
WD binary systems with various metallicities and IMeng & Yand (|2009|) presented the distribution of 
all the parameters for these close systems before the RLOF and at th e moment of SN la explo sion. 
Incorporating their results into the binary population synthesis code of iHurlev et al.l d2000l |2002|) . we 
obtain the distribution of the wind mass, M w j„d (see Fig.Q]) and the c olor excess via equ ations[T]and|6l 
The basic parameters for Monte Carlo simulations are same to that in lMeng et alj d2009l) when primor- 
dial binary samples are generated. Because the code is valid just for Z < 0.03, only five metallicities 
(i.e. Z = 0.03, 0.02, 0.01, 0.004 and 0.001) are examined here. 

The greatest uncertainty of our model is from Vwind- Here, we assume that Kvi n d = 10 km s _1 and 
we will discuss it in section |4] 
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Fig. 3 The distribution of color excesses of SNe la at maximum light for Z = 0.02 and 
various wind velocities (in km s" 1 ). The solid hi stogram is from observation, and the bar 
represents its maximum error dReindl et al.l 12005*1) . Via K-S test, the percentages show the 
confident level that the distributions of color excess from theory and observation are indistin- 
guishable. The black dashed line is the best fitted line, where the weights for the components 
of Z=0.01, 0.02 and 0.03 are 20%, 10% and 70%, respectively. 



3 RESULTS 

In Fig. |2] we show the distribution of the color excesses of SNe la at maximum light for v arious metal- 
licitie s. The distribution from observation is also shown by solid histogram in the figure dReindl et al.l 
120051) . We see in the figure that the SD model with optically thick wind can reproduce the distribution 
of color excesses of SNe la observed. The K-S test shows that the cases of Z = 0.02 and Z = 0.03 
have the highe st confident level th at the distributions from theory and observation are indistinguishable. 
The sample of lReindl et alj d2005l) includes the SNe la with various metallicities. H owever, it is difficul t 
to determine the m etallicity of the host galaxy of every SN la for a large distance dHamuv et al.ll2000l) . 
IWang et ail d2006l) collected the properties of several SNe la from previous papers and noticed that the 
mean metallicity of the host galaxies of the SNe la is [12 + log(0/H)] moan = 8.85 ±0.10 (private 
comm unication), which is consistent with solar metallicity ([12 + log(0/H)] so i ar ~ 8.8 dZaritskv et al 



1994). Then, our result is consistent with observation although the sample collected by I Wang et al 
d2006h is small. 

There are some SNe la w hose color e xcesses are very large (i .e. larger than 1. 0), such as SN199 9cl 
Jjhall2002t [Reindl et alj|2005l). SN 2003cg (lElias-Rosa et al.ll2006l) and SN2006X dWang et al.ll2008b . In 



the sample of Reindl et al.l ( 20051) . only one among 113 SNe la has a color excess larger than 1. In our 



simulation, the possibility of the high-color-excess SNe la is 2.0% for Z = 0.03, 0.4% for Z = 0.02 
and for Z < 0.02, consistent with that of the observations. This result might imply that the SNe la 
may have a color excess larger than 1 only when their host galaxies have a metallicity larger than 0.02. 
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Observationally, t he host galaxies of SN1999cl and SN2006X , i.e. NGC4501 and NGC4321, are both 
oversolar galaxy dCaputo et al.ll2000t iDors Jr & Copettill200r3) . Although there is no information about 
the metallicity of the h ost galaxy of SN200 3cg (NGC3169), NGC3169 is a Sa galaxy and probable has 
a oversolar metallicity dWillner et al.ll 19851) . 



4 DISCUSSION 



Our analytic model is so simple. We discuss various uncertainties about our model in this section. 



4.1 wind velocity 

The major uncertainty of our model is from the assumption that V w { n( i = 10 km s _1 . Although the 
consistency between the theoretical distribution of color excess and that from observation upholds 
this assumption, there i s no direct observational evidence to verify it. Many observational efforts were 
projected to find C SM dHamuv et al.ll2003l lAldering et aT1l2006t iPanagia et al1l2006t lOfek et alj|2007t 



Patat et al1l2.007 allbb . and only one observation obtained the constraint of the wind velocity ( Patat et aTl 
2007ah . The wind velocity is constrained to be smaller than 50 km s 1 (see Fig. 2 in lPatat et alj|20073h 



and then Patat et al. suggested that the progenitor of SN2006X should be a WD+RG system. However, 
we can not rule out the possibility that the observed CSM is fro m the opti cally thick wind since the 
companion has not been directly observed dHachisu et al. 2008). Recently, iBadenes et al.l d2007l) ex- 
plored the relationship between the SD models with optically thick wind for SNe la and the properties 
of the supernova remnants that evolve after the explosion. They found that the optically thick wind with 
velocity larger than 200 km s" 1 would excavate large low-density cavities around the progenitors. The 
large cavities are incompatible with the dynamics of the forward shock and the X-ray emission from the 
shocked ejecta in all the SNe la remnants as they examined. However, they also showed that if a wind 
velocity of 10 km s _1 is adopted, the properties of type la supernova remnants are well compatible with 
the prediction from the SD model with optically thick wind. Generally, the escape velocity from a white 
dwarf is in an order of m agnitude of 10 3 km s _1 , which is upheld by observations from rec urrent novae 
dWood & Locklevl2000l) . In the theoretically framework laid down bv lHachisu et alj dl996|). an optically 
thick w ind is formed in CO WD envelope with photospheric velocity of ~ 10 3 km s" 1 dHachisu et al.l 
Il999albh . where the expansion o f the photosphere is driven by helium flash in a helium shell on top of 
CO WD dKato&Hachisulll999t) . Then, in both theory and observation, it seems not to uphold an as- 
sumption of a low wind velocity at present. We check the influence of wind velocity on the distribution 
of color excess for the case of Z = 0.02 and the results are shown in Fig. [3] We see in the figure that 
when V > 100 km s~\ K-S test gives a low confident level that the distributions of color excess from 
theory and observation are indisting uishable, while when V < 50 km s _1 , K-S test shows an acceptable 
level. This result is similar to that in lBadenes et al.l d2007l) . 



4.2 delay time 

The secondary uncertainty of our model is from t^. Since optically thick wind may stop before SN la 
explosion, used in this paper overestimates the delayed time of the wind for some systems. However, 
this is not a serious pro blem since for most cases, SNe la occur during t he wind phase or after the wind 
phase for a short time dHan & Podsiadlowskil 120041; iMeng et ai1l2009l) . In addition, obtained from 
equation (10) only approximates the mean value for a given M\. This is also not a serious problem since 
we only check an average distribution of color excess. Generally, the peak of the distribution of color 
excesses moves to lower color excess with t^, and the peak value increases while the percentage of high 
color-excess SNe la decreases with t^. 
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4.3 metallicity 

The combination of different metallicities is checked and shown in Fig.[2]by the black dashed line since 
not all SNe la are from a unique metal licity. The we i ghts o f the components for Z= 0.01, 0.02 and 0.03 
are simply fitted from the samples of lHamuy et al.l ( 2000l) and IWang et al.l (120061) . We also show the 
best fitted line in Fig.[3]by the black dashed line. We can see from these lines that the combination of 
different metallicities does not significantly increase the confident level that the distributions for theory 
and observation are indistinguishable. Actually, for any combination of Z = 0.01, 0.02 and 0.03, the 
results will be acceptable, i.e. yielding similar K-S test probability, which is derived from similar K-S 
test probability for the distributions of Z = 0.01, 0.02 and 0.03. So, the basic results here still hold. 

4.4 interstellar extinction 

While we are not su re about the existen ce of the optically thick wind (it is a prediction by a model, 
not an observed fact, lHachisu et aill 19961) . we are absolutely sure about the existence of the interstellar 
extinction (IE). In the paper, we did not consider the influence of IE since it is difficult to separate it from 
CSM dust. However, we know that Type la SNe, as any other celestial object, suffer from interstellar 
reddening, arising in material that has nothing to do with the circumstellar environment of the exploding 
star. Let's consider the case of SN 2006X, whose color excess is certainly larger than E(B-V) =1, while is 
most likely arising in a cold molecular cloud which has nothing to do with the explosion site dWang et al.l 
120081) . So, for the case suffering interstellar reddening, we should remove the influence of IE on color 
excess. But at present, it is very difficult to do this. However, it is a clear effect that IE will lead to certain 
directions of where parameters such as Vwind and V c j cc should go, e.g. moving some very reddened 
observed SNe la to lowering reddened one would yield a higher Vvvind or y c j CC . For example, if we 
assume rather arbitrarily and simply that all observed SNe la suffer a extinction of E(B — V^)host = 0.1 
within its host galaxy, our model suggests Vwind = 20 km s _1 or = 20000 km s . Similarly, if 
E(B — y)host = 0.2, our model suggests Vwind = 50 km s _1 or V^j ec = 50000 km s _1 . So, the effect 
of IE could be counteracted by another uncertainties, V e j e c, an d tne basic results in this paper are still 
valid, (see the discussion about ejecta velocity below). 

If the color excess is from CSM as shown in this paper, one may argued that we would have seen 
the CSM material through radio and X-ray emission, ar ising in the shock produced by the fast moving 
SN ejecta crashing into the CSM dStockdale et alj[2~006l) . However. iPanagia et al.l (|2006j) showed a non- 
detection result for radio for 27 SNe la via the Very Large Array (VLA) observations, even including 
SN 2002ic, which indicated that mass-lose rate should be lower than ~ 3 x 10~ 8 M^yr . We take 
a typical case to chec k whether our mo del contradicts with observations. We set companion mass to 
2.1M Q (see Fig .10 in lMeng et al]|2009l) . which corresponds to a delay time of 2.5 x 10 6 yr. The mass 
of lost material is set to O.3M (see Fig.Q]). The mean mass-lose rate is 12 x 10~ 8 M Q yr _1 , which is 
higher but still comparable to that inferred fr om observations. A t tempts to detect radio emission from 
SN 2002ic with the VLA were unsuccessful dBerger et alj|2003t IStockdale etafll2003l) . SN 2002ic is 



the first case to show the signal of CSM flH amuy et al.ll2003l) and the amount of the CSM is 0.5-6M, 



dWang et al.ll200llChugai & Yungelsonll2004HUenishi et al.ll2004t iKotak et ai]|2004j) and its properties 



may be explained successful by the WD+MS mo del used here ( Han & Podsiadlowski 20061) . In addition, 
another r ecent twins of SN 2002ic (SN 2005gj. lA.ldering et aL 20061) has also not been detected using 



the VLA dSoderberg & Fraill2005i) . These non-detecti on results may indic ate that the mechanism that is 
successful used for SN Ib/c may not work for SNe la dPanagia et al.ll2006h . So, our model is at least not 
inconsistent with observation at present. 

4.5 ionization and ejecta velocity 

Generally, when a SN explodes, its radiation field is quite strong, especially at distance rjj < 10 16 
cm (for SN 2006X, rjj < 4 x 10 15 , |Patat et alj|2007al) . This radiation can easily ionize hydrogen and 
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Fig. 4 The distribution of color excesses of SNe la at maximum light for different ejecta 
velocity, where wind velocity is assumed to be 10 km s _1 . The ejecta velocity is unit in 10 4 
km s _1 . The black dashed line is the best fitted line, where the weights for the components of 
Z=0.01, 0.02 and 0.03 are 20%, 10% and 70%, respectively. 



evaporate dust up to quite large distances, causing the disappearance of reddening. In our simple model, 
this effect is not considered completely. We check the influence of the effect. 

We set ?'h = 10 16 cm and assume that all hydrogen atoms are ionized within the shell with radius 
of rn- The column density of ionized hydrogen is ?ih,o = / ri H ar~ 2 dr, where r\ = 1.73 x 10 15 cm (a 
ejecta velocity of 10 4 km s _1 and a rise time of 20 day are assumed). Then, the relative uncertainty from 
the ionization effect is = ^ H ' ~ r , l)rm ~ rH '°~ ri ~ 0.83 (r m > n). The uncertainty seems too 

>ih (r m -ri)r H ,o m,o 

large. However, this situation may be improved by taking a higher ejecta velocity and a smaller distance 
rn- A distance of rjj = 4 x 10 15 like SN 2006X may reduce the uncertainty to 0.57. 

The ejecta velocity used here ma y be smaller than the terminal velocity of SN la which can be as 
high as 3 x ^kms" 1 dWandl2006l) . This value may reduce the uncertainty from the ionization effect 
to ~ 0.48, which means that no more than a half of CSM are evaporated up to large distances and have 
no influence on the color excess of SNe la. Then, the ionization effect at most decreases the confident 
level that theory and observation are indistinguishable from 74% to 65%. 

In Fig. |U we show the influence of ejecta velocity on the distribution of the color excess of SNe 
la. The influence is similar to that of wind velocity. From the figure, we can see that the peak of the 
distribution increases with the ejecta velocity, and when V c j cc < 5 x 10 4 km s _1 , the distributions of 
color excess from theory and observation are indistinguishable on the confident level of higher than 
50%. Fortunately, the effect of ionization and a higher ejecta velocity is same to that of interstellar, and 
their influence on the decrease of the confident level may counteract each other. Then, we can say that 
the distribution of color excess between theory and observation are indistinguishable on the level of 
higher than 59%, even when y c j CC = 3 x 10 4 km s _1 . 
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For the discussion above, our model may get a meaning result although our analytic model is so 
simple. 

5 SUMMARY AND CONCLUSIONS 

In summary, if a wind velocity of ~ 10 km s _1 is adopted, the SD model with optically thick wind may 
reproduce the distribution of color excess of SNe la obtained from observation, which might support the 
SD model for SNe la. However, if a wi nd velocity larger than 100 km s _1 is adopted, the reproduction is 
bad. Our results are similar to those of lBadenes et al.l d2007l) . Then, it should be encouraged to perform 
more detailed observations about the velocity of CSM around progenitors of SNe la. 
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